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In the work presented here different models employed 
to describe ferromagnetic resonance (FMR) and ferromagnetic 
relaxation (hereafter referred to as relaxation) have been 
briefly discussed. The " s-d exchange" interaction model 
between ferromagnetic 3^ electrons and conduction s electrons 
in a ferromagnetic transition metal (Ni, Fe) (hereafter referred 
to as ferromagnetic metal) has been studied m detail to explain 
relaxation m these metals. The coupling of spin waves (magnons) 
with lattice vibrations (phonons) resulting from magnetostriction 
has been considered and a two magnon - one phonon process has 
been studied. The renormalized spin wave and phonon energies 
have been obtained. The method of double -time temperature 
dependent Green functions has been used to investigate these 


effects 
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Chapter I gives a brief account of the present state 
of our knowledge of FMR and relaxation, both m ferromagnetic 
insulators and metals. FMR and relaxation m ferromagnetic 
insulators is well understood m terms of various magnon-magnon 
and magnon-phoncn processes However, the lmewidths and 
relaxation times m ferromagnetic metals are not explained 
on the basis of these processes In a ferromagnetic metal 
the lmewidth is believed to be due to (i) magnetic relaxation 
damping - the microscopic origin of which is not exactly known 
This is described macroscopically by including Landau-Lif shitz 
or Bloch-31oembergen type of damping terms m the equation of 
motion for the magnetization of the system ( 11 ) exchange 
conductivity broadening - which arises because of the fact 
that the exchange stiffness opposes the tendency of the 
dynamic magnetization to vary within a small skin depth. The 
lmewidths and relaxation times in ferromagnetic metals are 
not explained on the basis of these two mechanisms even 

In Chapter II Kittel's phenomenological theory of FMR 
has been briefly discussed The expression for the FMR 
frequency has also been derived quantum mechanically as 
obtained by Van Vleck The effect of the shape of the 
specimen, magnetic anisotropy and domain structure on the 
FMR frequency has been qualitatively discussed The 
phenomenological approach to damping has also been briefly 
described. The phenomenological theory satisfactorily 
explains the shape of the resonance curve but is inadequate 
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to explain lmewidths and tnexr temperature dependence It xs 
not surprizing as the microscopic nature of the relaxation 
mechanism is not exactly known 

Chapter III presents a brief introduction to the Green 
function formalism that has been used m this and the subsequent 
chapter The method is applied to the Heisenberg model of a 
ferromagnet with demagnetizing energy as first considered by 
Tyablikov . It is noticed that m the approximation (random 
phase approximation) considered here the Green functions have 
poles on the real axis and consequently there is no damping. The 
FMR frequency formula obtained is applicable at all temperatures 

The coupling of the magnons with phonons has been 
considered and a two magnon-one phonon process has been studied, 
The renormalized magnon and phonon energies are obtained 

In Chapter IV the spin waves in an itinerant electron 
model have been discussed. The interaction of the ferromagnetio 
3d electrons interacting with conduction s electrons by means 
of "s-d exchange" interaction has been considered The express- 
ions for the dynamic magnetic susceptibility and for the line- 
width due to "s-d exchange" interaction have been obtained. It 
has been suggested that a bunch of magnons participate in the 
relaxation process resulting in the conduction electron spin 
flip. 

In the concluding Chapter V, the results of the present 
investigation are summarized. The limitations of the theory 


are also discussed* 
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CHAPTER - I 
INTRODUCTION 

Ferromagnetic resonance (FMR) is a powerful experimental 
tecnnique for studying the modes of excitations of a ferro- 
magnet and the lifetimes of modes arising out of their coupling 
with other excitations The basic principle of the experiment 
is as follows 

A ferromagnet is placed m a static magnetic field and 
an oscillating radiation field is applied perpendicularly to 
the constant field*. As the frequency of the radiation field 
is varied, resonance occurs due to absorption of energy from 
the radiation field at a particular frequency depending on the 
value of the constant field. The ferromagnetic resonance of 
this type was first reported by Griffiths 1 m 1946. Jt differs 
from other kinds of magnetic resonance such as electron para- 
magnetic resonance in that here the collective effects are 
involved. All the spins are parallel and collectively precess 
about the static field in phase It is also called uniform- 
mode ferromagnetic resonance . In the framework of the spin 
wave theory the uniform precessional mode corresponds to 
excitation of spin waves of wave vector k = o . (The excitation 
of spin waves of k £ o requires rather special conditions. 
These non-uniform (k /£ o) modes are excited in spin wave 

*There is a different type of experiment m which r.f. field 
is applied parallel to the static field. This is called a 
parallel pumping experiment. 
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resonance experiments These experiments were suggested by Kitte. 

2 

and first observed by Seavey and Tannenwald } . 

In the present work we would restrict our discussion 
to the uniform mode ferromagnetic resonance. The occurrence 
of this resonance is due to the transition between the state 
of the system with total spin NS (where N is the number of 
lattice points and S is the atomic spin) aligned along the 
magnetic field direction and the state with all the spins 
reversed, i.e., with total spin - NS. The transition is 
caused by the exciting radiation field. 

The first experiments on FMR done on ferromagnetic 

metal (Ni, Fe, Co) thin films by Griffiths^ indicated that 

the absorption curve had a resonance behaviour with a maximum 

for a certain value of the magnetic field which depended on 

the frequency. The interpretation which he gave suggested 

that the resonant absorption by the magnetic dipoles (which 

might, of course, be multiples of where -q 1S 

Bohr magneton) was taking place leading to a loss of energy 

by the radiation field. The resonant frequencies, he observed, 

were greater than the calculated Larmor frequencies for 

electron spin by factors of about two to six. Griffiths 

unsuccessfully tried to explain this anomaly by introducing 

Lorentz field, an approach which was not justified. A correct 

3 

theory was later given by Kittel^ m 19^7 which led to the 
values of the resonant frequencies in close agreement with 
experimental determination. Kittel showed that it was 
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important to consider the dynamic coupling caused by the 
demagnetization field normal to the surface of the specimen. 

The principal result of ms theory was that the appropriate 
Larmor frequencies should be calculated as for a field 

i 

strength (BB^)^ rather than for B, the static magnetic field, 
where B^ is the magnetic induction 

The width of the resonance line arises out of the 
interaction of the uniform mode with other modes (such as 
magnons, phonons and electrons) m the system The origin 
of the width of the ferromagnetic resonance lines is not 
clearly understood specially m metals 

We give below a brief account of some experiments on 
FMR which mainly deal with the measurements on lmewidths 
and relaxation times and their interpretation both in insulators 
and metals 
Insulators 

In non-conducting f errimagnetic single crystals FMR 
(kjv o) lmewidth can be very narrow. LeCraw and Spencer^ 
who measured lmewidths m pure yttrium iron garnets (YIG) 
found that the lmewidth can be as small as 0 2 oersteds at 
3*3 KMC/sec 

The dominant relaxation mechanism in these materials 
appears to be magnon-magnon interactions. For a detailed 
discussion of the various multiple magnon processes we refer 
to Haas and Callen^and to Sparks^. 
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In any resonance experiment one might expect some 

broadening of the resonance line due to specimen inhomogeneities 

7 

or surface imperfections LeCraw et al showed that the 
ferromagnetic resonance lmewidths of yttrium iron garnet 
spheres depended on the size of the grit used for polishing 
the specimens and hence to the surface inhomogeneities Sparks 
et al have proposed that the two magnon scattering of k p o 
magnons plays an important role m the relaxation of imperfectly 
polished YIG. 

Q 

Spencer and LeCraw^ have measured the longitudinal 
relaxation time ^ m YIG as a function of temperature. 

They find that at 9»3^ KMC/s between 2° and JO°K 

\T\ io 5 T sec"* 1 

w 1 

which is in good agreement with the value calculated by Sparks 
and Kittel 1 ^ 

Various types of volume inhomogeneities can also affect 
1 1 

the line width 

In principle, the magnon-phonon interactions should 
also contribute to the line width and a direct process would 
occur when the phonon and magnon frequencies and wavelengths 
are matched In the "cross -over" region 'where magnon and 
phonon frequencies and wavelengths are equal) the magnons 
experience large magnon-phonon interaction and this results 
in a marked increase jn the of these magnons wr.ere 4f^k 
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is the width of the spin wave resonance. An increase m /\ 

will be reflected m an increase m the critical r.f field 

strength necessary to initiate unstable growth of these spin 

waves Turner ^has found that m YIG the critical r.f. field 

m a parallel pumping experiment (B i 1 b x ) does indeed increase 

sharply for certain values of the d c. field, that is, for 

13 14 

certain values of k. Akhiezer and, Kittel and Abrahams 
have theoretically investigated phonon-magnon processes m 
detail. 

Yttrium iron garnets containing rare earth impurities 
have also been experimentally studied. The FMR lmewidths 
of rare earth garnets and of YIG containing rare earth 

■i e 4 

impurities differ markedly 0 ’ from those of the pure YIG 

m that the lmewidths of the rare earth garnets are very 

large In addition, the lmewidth has a very strong temperature 

dependence . The lmewidths have been explained on the basis 

of magnon-phonon scattering induced by rapidly and slowly 

17 1 8 

relaxing rare earth ions * * 

Relaxation processes involving transfer of valence 

electrons between ions of variable valency such as ferric 

5 

and ferrous ions have also been studied m garnets. 

One can fairly claim that the lmewidths m ferromagnetic 
insulators and ferrites are well understood m terms of 
magnon-magnon and magnon-phonon processes, which are induced 
by dipolar interactions and magnetostriction. 
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In the case of metals, however, the situation is not 

very clear. The electromagnetic waves enter a narrow surface 

layer of the metal called the skin depth and consequently the 

direction of magnetization vector M is not uniform and the 

effective exchange field is not parallel to M and a torque 

acts on M which modifies the equation of motion of M and 

the nature of the resonance. This modified resonance is called 

i q 

11 spin wave resonance" following Rado and Weertman . As stated 

earlier, these spin wave experiments were first suggested by 

20 2 
Kittel and rfere first observed by Seavey and Tannenwald 

In very carefully prepared specimens of Co, Ni and Fe, 

the narrowest spin wave lmewidths that have been observed by 
21 22 

Rodbell , Ihillips and Rosenberg are of the order of 30 - 50 
oersteds. 

Tne lmewidth measurement s on perfect metal single 
crystals have been made by several other author s 2 ^ ’ 23> ^ 

The relaxation times obtained from their measurements cannot 
be explained on the basis of spxn-spin^’^ and direct spin 

•i o piy pQ 

lattice relaxation ’ ’ . They, however, show that the 

relaxation of the spin waves excited at FMR contributes to 
the lmewidth m addition to the exchange broadening caused 

PQ 

by mhomogeneity of magnetization in the skin depth Even 
if surface pinning of spins is included, this result remains 


unaffected 


30 


More recent and detailed measurements of linewidths 
on pure iron are by Bhagat et al" 31 ’ 32 . The widths in iron 
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31 

are found to vary as the square root of the frequency and 

are due mainly to exchange conductivity mechanism. They were 

also led to conclude on the basis of their observations that 

for iron the relaxation frequency /\is essentially temperature 

independent from 4°K to 1000°K 

In the case of Ni , more recently lmewidth measurements 

have been reported by ilodbell-^, sa+amon^, and Bhagat et al^ 2 ’^. 

The lmewidths m Ni are enormously large (of the order of a 

few thousand oersteds) The results obtained by these authors 

are not all m agreement The widths for Ni can be accounted 

for by relaxation damping described by a frequency independent 

Landau -Lif shit z coefficient )\(= 2.3 • 1 0 8 sec’ 1 ) plus a 

smaller contribution from the exchange conductivity mechanism^ 1 

Salamorr?^ finds an effective (relaxation frequency) 

decreasing with increasing temperature. This result is xn 

qualitative agreement with Bhagat et al . Salomon also 

finds an increase in g value as he goes from below T to 

c 

above T c> In the case of Bhagat et al, the increase m g 
is rather gradual. The reasons for this discrepancy are 
not clear . 

In order tc understand the mechanism of FMR, it is 
necessary to consider various models that have been used to 
describe it The important models are the following 
Spin Wave hode! 

This model assumes that the system is very near to 
its ground state. It is, therefore, particularly valid at 
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low temperatures when the departure of the ferrcmagnet from its 
ground state (when all its spins are aligned parallel to each 
other } is very slight. Van Kranendonk et al J used t.iis model 
to obtain tne frequency of the uniform ferromagnetic resonance 

o 

fhe explicit structure -5 of the Hamiltonian used is the 
following 



where Z 

D _ , r denrtes the vo-ctor joining th~ atoms 

J r J / ° x and 3, 

ij /i B is the Bohr magneton, B is the 

constant magnetic field, J denotes the isotropic exchange 

J.U 

integral, denotes the atomic spin at the 1 lattice site. 

The various terms m equation (l) represent respectively 

the Zeeman energy, Heisenberg's exchange energy and the dipolar 

energy. 


The interaction of the radiation field with the spin 
system assuming that the radiation field B x is uniform 
tnrougnout the sample (i.e., the skin depth is longer compared 
to the dimensions of the specimen under investigation) has the 
form given by 


H int = $ x 


h «* £ s ± 


( 2 ) 


where the sum is over all the lattice sites of the crystal. 
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The radiation field thus couples only with the transverse 
components of the magnetization 

The interaction (2) causes transitions between different 
energy levels of the ferromagnetic electron states and excites 
spin waves of wave vector k = o (called uniform precession 
mode) as can be seen by replacing by in tbe 

harmonic oscillator approximation, q being the coordinate 
operator . Thus 

*U* %. * SyU^NS)* % 

1 TT-l 

where q Q = — > q is the completely symmetrical normal 

coordinate. Thus the radiation field excites only the uniform 
processional mode. 

36 

Van Kranendonk et al using harmonic oscillator 
approximation obtained the normal modes of the system The 
spin waves w.uch are the elementary excitations are the normal 
modes of the spin system The energy eigenstates o± the system 
are characterized by n^, the number of excited spin waves 
wit a. wave number k. Since the k = o spin waves are excited 
at FWR, the frequency of this mode is given by 

lOo * S/' b |_[b + ( 3 * (V"* )**■]]* (4) 

This result is identical with Kittel* s expression except 
that the magnetization M is replaced by the saturation 
magnetization value to . This is expected as the spin wave 
model gives the low temperature value of the FtoK. frequency. 
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The model Hamiltonian equation (l) is based on the 
assumption that the electrons are localized on the atoms at 
the lattice sites It is particularly suited to the case of 
ferromagnetic insulators where the localized spin picture is 
applicable The localized spin model may also be used to treat 
rare earth metals because m them the 4f shell is well within 
the atom and the electrons can be regarded as localized The 
situation m ferromagnetic transition metals is difxerent 
Here the 3d shell is on the outside of the atom and the 3d 
electron wavefunctions of the neighbouring ions overlap and 
thus form a narrow band of states This means that the 3d 
electrons move from one atom tc the other and, therefore, the 
localized electron concept is not applicable This is called 
an itinerant (band) electron model"^ - ^. This is based on the 
competition between the kinetic energy of the electrons m a 
band and the electron exchange m the Hartree-Fock approximation 
This exchange may be thought of as arising from the intra-atomic 

1 4c , 

rather than inter-atomic exchange. Ilajagopal et al nave 
studied the magnetic properties of itinerant electrons m 
detail. We shall discuss spin wave excitations m this model 
in Chapter IV 

Various modifications of these itinerant and localized 
models have also been used to treat real solids One such 
model is the so called "s-d exchange" model. This has been 
used to investigate lmewidths and relaxation times in 
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ferromagnetic metals It was first considered by Zener 
42 

Kasuya showed that the exchange interaction between the 3d. 
electrons is mediated by s electrons and leads to an effective 
Hamiltonian for d electrons The Hamiltonian for the magnetic 
3d spins exchanged coupled to s band electrons can be written 


m the form 


43,44 


H = H + H , 
o s-d 


H = 
0 

Xc c + 

„ q q 

c 

q 



q 


m /£ n 

and 

H . = 

s-d 

- 2 2! 

J s-d (V 



mn 


where £ 

q 

rO 

3 O — 

q 

3 J 3 -d( qq ) 


J (m - *n) * 2 Ag B £ S s= 

/ n n 


c and c 
q q 


3 S 
m n 


where £ q = £ q - s J g _ d (qq) + 2^^ B (5) 

*2 2 

C° s is the renormalized s electron energy 

q 2m 

c + and c denote creation and annihilation operators 

q q 

for s-band electrons. J Cm -~n) and J . denote the d-d 

and s-d exchange integrals. In this model also the 3d electrons 

have been assumed to be localized which is not justified as 

has been pointed out earlier. 

By writing the Hamiltonian ( 5 ) m second quantized form 

43 44 

and. also using spin wave approximation, it can be shown ’ 
that the "s-d exchange" interaction brings about the coupling 
of the spin waves and conduction electrons which are the 
elementary excitations of the solid m this model The spin 
waves can relax through their interactions with conduction 
electrons if they (spin waves' have sufficient energy (of the 
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order of J g required to flip the conduction electron spin) 
This model has been examined m Chapter XV of the present work 
m detail. 

The spin waves can also couple with lattice vibrations 
(phonons). This magnon-pbonon coupling can be brought about 

1+5 46 

by the magnetic dipolar interaction and/or magnetostriction 
A two magnon - one pnonon process has been studied m 
Chapter III. 
jr-henomenological Model 

Je know that the equation of motion of the magnetization 
vi can be written as follows 

# - n* * (6 > 


where is the spectroscopic splitting factor. This 

is the undamped equation of motion for f?. To explain 
resonance lmewidths it is very essential to include damping 
into the equation of motion. The first phenomenological 

47 48 

equation was suggested by Landau and Lxf shitz and Gilbert 


and reads as follows 


49 


dM 

dt 


= /[ M x B + y ^Anisotropy torque~J + M 

(7) 




— ^ d 

where and A are respectively saturation magnetization 

and exchange stiffness, y) is the relaxation parameter. The 
terms xn (7) respectively denote the torque due to external 
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(static) field, anisotropy torque, the exchange field torque 
and the spin lattice torque By working with a macroscopic 
dynamical equation one can obtain information about the 
important parameters. 

Due to skin effect we have two important regimes m 
metals - normal and anomalous. The lmewidth m the normal 

conductivity regime (1 , 1 is the mean free path of the 

r 31 

electron, ^ is the skin depth) can be roughly given by ( 31 ) 


1 ,45 7 ) 


C A t 6d 2 


where the first term is contributed by the magnetic relaxation 
damping and the latter by the exchange conductivity mechanism. 
It is evident from (8) that the lmewidth will be dominated by 
exchange effects at low frequencies and by magnetic relaxation 
damping at high frequencies^ 

Anomalous Skin Effect Region 

In the low temperature region where the average mean 
free path 1 becomes larger than the skin depth (i e, , 
the skin effect becomes anomalous In this regime the 


effectiveness of the electrons m absorbing energy from the 
radiation field is reduced, or stated differently the 
conductivity becomes k (wave number of the spin wave 
excited) dependent and decreases for increasing k. 



In the extreme anomalous skin effect region (l^<5) the 


conductivity becomes proportional to k. Since the value of k 

at resonance is large ({v_10^) this means that the effect of 

anomalous conductivity will be to reduce the exchange effects, 

as compared tc the case wnen the conductivity is taken equal 

4 3 

to its k = o value (normal region) Vonsovski et al and 
Turov^ find that the lmewidth becomes independent of 
temperature in tms region 

4.3 

Vonsovskal et al have shown that the spin wave damping 
can take place only when it has sufficient energy (of the order 
of J q ) required to flip the conduction electron spin. The 

K 

energy of a typical spin wave (k Or ) is of the order of 

—6 ^ 

10 eV. This is several orders of magnitude smaller than 

- 1 42 

J g d ~ 10 eV found by Kasuya . Therefore, the spin wave 
damping will not be favoured for low energy spin waves due to 
the large exchange gap between bands of s-electrons with 
different spin orientations and the energy conservation. 

44 

Recently , however, it has been pointed out that this in 

itself is not such a serious difficulty because m principle 

non-resonant aosorption of energy to non- stationary states 

of s-spins is possible at very low magnon energies Kittel 
<5 1 

and Mitchell considered the effect of s-d exchange interaction 

on the lmewidth. They assumed that the conduction electron 

is relaxed to the lattice independently by a rapid mechanism. 

Turov'* 0 has shown that the lmewidth depends very strongly 

on J , ( £>£ The value of J . is not precisely 

s-d — s-d' s-d 
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known However, to find agreement with experiments a value 

_2 

of J g of the order of 10 eV wnl be required This 

value is not entirely unreasonable but is an order of magnitude 

42 

smaller than that found from some other measurements 

Secondly to fit the frequency dependence of the lmewidth a 

frequency dependent J g is required which is highly 

43 50 

unlikely. According to Vonsovskn. et al and Turov the 

lmewidth should decrease, (m the normal conducting regime) 

with decreasing temperature m direct proportion to electrical 

resistivity This does not agree with the observed behaviour 

4q 50 

m any temperature range . Turov also considered the 
contribution made to the resonant lmewidth by the interaction 
of the s-spins with the magnetic field due to the spin waves 
and showed that these processes made a much smaller contribution 

44 

to the lmewidth Recently Heinrich et al has investigated 

the influence of s-d exchange on relaxation of magnons m 

metals (Ni, Fe). By phenomenologically taking into account 

the finite value of the s-spm lattice relaxation time they 

seem to be able to account for their observations on temperature 

24 

dependence of the relaxation frequency /) m iron . However, 

31 52 53 

the measurements of Bhagat et al ’ ’do not show any such 

dependence and, therefore, the applicability of their theory 

to iron is far from being unequivocal. They, however, admit 

that their theory is not likely to explain the observed 

54 

lmewidths m nickel. Very recently L®ichner et al have 
claimed to have been able to explain the observed linowidths 
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55 

in KgCuCl^ 2HgO using Van Vleck‘ s moment expressions 
Unfortunately, it is not possible to explain the Imewidtiis 
m metals (specially Ni) on the basis of existing theories. 
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CHAPTER - II 

I- HENQ^JS NGLOGI CAL THEORY OF FERROMAGNETIC RESONANCE 

The results of measurements on resonance frequency and 
lmewidth m FMR can be interpreted xn terms of a macroscopic 
dynamic equation called, the phenomenological equation^’^ and 
information about important relaxation processes m a ferro- 
magnet can be obtained from such an equation The phenomeno- 
logical theory of FtoR is thus a macroscopic theory. 


2*1 Effective Intern a l Field m a Ferromagnet 

The equation of motion for the magnetication vector it 
of a ferromagnet placed m a static magnetic field B in z- 


dir action is of the form 


= Y f_M x H "j 


( 2 . 1 ) 


/ ge 

= 2mc ' 1S ^'^ le spectroscopic splitting factor and 
is the effective internal magnetic field inside the 
ferromagnet. It includes besides the external field all 
internal fields arising out of magnetic anisotropy energy, 
demagnetizing dipolar energy and magnetostriction (state of 
stress) Lorentz and "Weiss fields are omitted as they are 

' — ^ y 

parallel to Ii and hence their cross product with to will be 
identically zero Equation (2.1) is the fundamental equation. 
We will use this to obtain the expression for the resonant 
frequency. 
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2 2 General Expression for the FMR Frequency 

3 

The FMR frequency was first obtained by Kittel using 
the phenomenological theory The effective values of the 
magnetic field components inside the ferromagnet, ignoring 
anisotropy effects, are 


H 


B - N M 

XI 


X XX 

li 

— 

-N i'. 

yi 


y y 

H 

jj. 

B - N M 

Zl 


z z 

: » N 

and 

N are called the 

x ’ y 


z 

M , 

M 

and M are the x, 

x’ 

y 

z 


(2.2a) 

(2.2b) 

(2.2c) 


factors M x , My and M z are the x, y and z components 
of the magnetization. The radiation field B x is along x- 
direction 

Since M and M are small compared with as the 

x y 

specimen is magnetized to saturation m the z-direction, we 
can safely neglect products such as MM in the equations 

* y 


(2.1) which become 


dh. 


dt 

dM 


~ = /[B + (N y - N z )M a 3 M y 


dt 

dM 

j 

dt 


z - /t M z B x -< N X - N z> M * fI z - V 3 


0. 


(2.3a) 

( 2 . 3 b) 

(2.3c) 


o 1 to t 

These equations ( 2 * 3 ) can be solved by setting M x = M x e , 

M x Cj * i~) 

M 5 M° and the susceptibility can be found 

Y Y x a x 

to be "Y = ° " f r — 5— where 7(~ is static 

/v, x . / uj . « 0 

1 ~ ^ t#o) 
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susceptibility and is the ferromagnetic resonance 

frequency for an ellipsoidal sample given by 

t-o ■ k_i B + ( n x - n z )‘ 4 > \ b + ( v K * ) M iJ T (z - 5) 

Here u has been replaced by 1 , 

2 . 3 Quantum Mechanical Derivation of FMR Frequency 

The Fj>© frequency can also be derived quantum 

36 kk eg 

mechanically ' ’ considering suitable models of the 

ferromagnet Cne such model is the Heisenberg model with 

magnetic dipolar interaction first considered by Van Vleclt . 

This model assumes that (a) the electrons are localized on the 

atoms (bj only one electron per atom is present (c) the electrons 

are exchange coupled through an interaction of the type 

-J _ *S* , whdre J is an exchange integral, S. and 
J * J 1 J X 

S are the spin operators corresponding to the electrons 
J 

on tne atoms at the lattice sites 1 and j Apart from 
this exchange interaction, there is also a magnetic dipole- 
dipole interaction This dipolar interaction depends not 
only on the orientation of the spins with each other but also 
on their orientations with the line joining them. 

The model Hamiltonian is of the form 


fi = g 


fr ? % 


o Z 

°J + 


S 7a 'koVf, *L 


k 7 j 


^ 1 k 


Jk 


( 2 . 6 ) 
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rxere the terms on the right hand side of equation (2.6) 

denote respectively the Zeeman energy (associated with the 

static magnetic field m the negative z -direction', , exchange 

energy and the magnetic dipolar energy The coefficients 

D, have the value 
kj 


D. 


kj 


2 , 2 
e M b 


jk 


if the dipole-dipole interaction is of true magnetic origin 
The values of D, larger m magnitude can also occur as a 
result of spin-orbit coupling Such dipolar interaction is 
called pseudo -dipolar interaction and is of snort range 
character. i/e do not consider such pseudo -dipolar coupling 
here 


Je rewrite the Hamiltonian (2.6) as 


n = g 


^ 3 IT 0 


> 

J 


k>j 


V t 


■=- * Z A kj S k S 3 -3 Z. "kJ ‘jk v-k- 0 k/'=j-*Jk' 




(*•« ) 


where 




= J kj + D kj 


The equation of motion method is now used to obtain the 
resonance frequency. 

The equation of motion for the atomic spin S can be 
written as follows 


dS . 


dt 


i 

h 


Ch, sl 

J 


(2.7) 


where j^H, , denotes the commutator of with the 
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Hamiltonian (2.6) of the spin system. 

Using the well known commutation relations 


r „y — i 

L S j> S k J 


O z 

i a , s 

jk j 


etc 


equation (2.7) yields for the three components of the total 
spin 5 


H 


dS 

dt 


= D S + 3 D kj L^jk b j + ^ jk 

A 73 


+ 3 


Jk °3 j 


Ajk S k + ^jk S k^ 


(2.8a) 


dS J 

dt 


g^ 3 B S* + 3 £ D kj ( ^ Jk 3j + f) 3k ♦ 4 S*) 

' K ^ J 


( °(]k 3 k “ ^ik S k' 


Jk 


Jk j 

(2.8b) 


and 


a 


dS 

dt 


3 JZ S I + >^k s 3> 

( " ^nk S k + /^jk S k^ 


(2.8c) 


Here ^3^ and denote the direction cosines of the 

vector r , joining the atomic spins at the lattice sites i 

JK 

and j. S = & denotes the total spin of tne system. 

3 3 

It has been assumed here that the skin depth is large 
compared to the dimensions of the specimen so that the radiation 
field is uniform throughout the sample. The magnetic operator 


S'* which is responsible for the absorption can be written as 


,x 

5 

-X 
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The exchange energy co mutes with any component of the 
total spin and that is why it does not appear m eqn (2.8). 

In the above equations (2.8) one can neglect products 
such as S.^ on the basis that in a ferromagnetic medium, 

the specimen will be practicallv i agnetized to saturation 
along the field direction and the matrix elements of S and 

yr 2J 

S will be smaller than those cf S . In the products such 

as S* S, Z , we replace S Z by its mean value — — — rr where N 
J k 3 g 

is the number of atoms in the crystal 

Using these approximations and assuming that the atoms 
in the crystal are symmetrically situated (so that the terms 
such as 2 oC jk ^ etc vanish) , we obtain 


dT ■ fa B s y + (VV*/^ » s 1 


* fa 3 S X - (N z -N x )g ^ H 


( 2.8 ) 


rv o 


where the shape dependent demagnetizing factors are defined as 
follows 

2 

s = -L- z - ~ 3 + 4ZT 

* N k r\ 3 

and the corresponding expressions for N and N 

y z 

j£qn, (2.8 ) can be solved by setting 


t,x _x itdt 

b = b e 
o 


, _y 0 y lWt 
and S = S e , 

o 


to yield the Kit tel resonance frequency (2.5) 
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It should, however, be pointed out here that the present 

model is particularly suited to the case of ferromagnetic 

insulators suca as garnets, since the localized electron spin 

concept is applicable m their cases Tms model was employed 
55 

by Van Vleck to explain ferromagnetic resonance absorption 
* 

a.id lmewidth 

2.4 Effect of the Shape of the Specimen on FMR Frequency 

As we know that In a ferromagnet there are strong 

57 

internal fields arising out of magnetic anisotropy caused 
by the magnetic interactions between electrons, the effective 
internal field will depend on the symmetry of the crystal, 

the shape of the sample, the magnetization and its direction 
and on the direction of the external static field B m 
relation to the crystallographic axes From the expression 
(2.5) it is clear that the resonance frequency depends on the 
shape of the sample as ;he demagnetizing factors enter into 
its expression. w?e consider some special cases of eqn, (2.5). 

(a) I- lane For an infinite plane sample in the x-z direction 

N x = N z = 0, N y = 4 IT 

and CO Q = / [b (B + 4-rr M z >3 * = /(BB 1 )^ 

where B^ denotes the magnetic induction. 7e , therefore, 
notice that is equal to the Larmor frequency for a 

fictitious field (B B^)^. 

(b) Sphere , For a spherical sample all the demagnetizing 
factors will be equal, i.e., 

The linewidth is defined as the width of the resonance curve 
at the half power absorption. 
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N x = S y = N z = 4 A /3 

and Wo = 

i e., for a spherical sample is equal to the Larmor 

frequency. 

2.5 Magnetic Anisotropy and its Influence on PMR Frequency 

Ferromagnetic crystals have "easy" and "hard" directions 

of magnetization, i.e., it is easy to magnetize a crystal m a 

direction (called easy direction) than m the other (called 

hard direction) The excess amount of energy needed to magnetize 

a crystal in a "hard" direction than m an "easy" direction is 

57 

called the magnetic anisotropy energy . 

In a crystal tnere is an intense, inhomogeneous and 

asymmetric electric field (called the crystal field) produced 

58 

by the constituent ions As a result of the crystal field 

and the spin-orbit coupling, the electron orbits are distorted 
along particular directions and the magnetic fields associated 
witn these orbits have directional properties and the spin- 
magnetic moment orients itself along the resultant of the 
fields produced by the orbits and thus anisotropy originates. 

The spin-orbit coupling is a consequence of the partial 
quenching of the orbital angular momentum by the crystal field, 
tfhat happens m an asymmetric electric field as exists in a 
crystal is that the orbital angular momentum components are 
no longer constant (in a symmetric field one component usually 
taken as L z and the square of the total orbital angular 
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2 

momentum L are constants of the motion) In a crystal L z 

will no longer be a constant of the motion, although to a 

2 

good approximation L may be a constant of the motion irfhen 

i/ z averages to zero, the orbital angular momentum is said to 

be completely quenched However, the quenching is usually 

not complete ' 3 ’ and, therefore, the spin interacts with the 

orbital motion wmch in turn interacts with the crystal lattice 

by means of the electrostatic field. Thus the spin (or the 

magnetization) "sees" the crystal lattice and its geometry 

Tms is why the anisotropy depends on the crystal symmetry. 

Owing to the spin-orbit coupling, the spin vectors 

responsiDle for ferromagnetism feel slightly anisotropic 

electrostatic forces which connect the orbital angular momenta 

of different atoms and whose bonding effect is called orbital 

59 

valence Consequently, anisotropic dipole-dipole and 

anisotropic quadrupolar couplings (usually called pseudo - 
dipolar and pseudo-quadrupolar interactions respectively) 
are produced. These interactions are used to explain magnetic 
anisotropy m f erromagnets . Van Vleck was the first to 
introduce pseudo -dipolar interaction to explain anisotropy 
of single cubic metal crystals. As will be mentioned later, 
Van Vleck also considered a microscopic model of anisotropy 
with a quadrupolar coupling between the spins to explain 
resonance lmewidth and frequency shift. 

For cubic crystals anisotropy energy can be written 


m the form 
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fi an = M°< 2 P 2 -P 2 / 2 +1U Z ) + ^ 2 (i*/ Z (2.9) 

where and Kg are the first and second order anisotropy 

constants respectively and j/ , j^> and denote the direction 

cosines of the magnetization referred to the crystal axes The 
constants K^ and Kg can be determined from experiment 
for iron at room temperature 

K 1 = 4.2 x 10- 5 erg/cm-*, Kg = 1.5 x 10^ erg/ cm 3 (2.10) 

For crystals with a single ^referred axis, such as 
cobalt, the anisotropy energy may be of the form 

£ an = K 1 Sin2 ^ + K 2 Sxn4 0 

where 0 is the angle between the magnetization and the easy 
axis, higher order terms are usually neglected For Cobalt 
at room temperature 

K 1 = 4,1 x 10^ erg/cm^. Kg = 1.0 x 10^ erg/cnf* 

( 2 . 11 ) 

The anisotropy constants depend strongly on temperature 

The anisotropy energy influences FMR frequency through 
its dependence on the angles formed by the static field with the 
crystal axes At a fixed radiation field frequency the value 
of the constant magnetic field at which the resonance occurs, 
will be lowest near the "easy" direction and highest near 
the "hard" direction. This angular dependence has been 
confirmed by the experiments of Kip and Arnold^ 0 , 
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3 

Kittel considered the effect of the macroscopic 

anisotropy energy of the form (2.9) and found that its effect 

can be included by the proper alteration of the demagnetizing 

5 5 

factors used m (2.5) • Van Vleckr^ , however, considered the 
microscopic model of anisotropy with a quadru^olar coupling 
between the spins of the type 


* 7 , 


K 


-4 ,2 ,2 


jk jk ' j * jk ; 


( 2 . 12 ) 


it should be noticed that for identical experimental situations 
both Kittel and Van Vleck get almost identical results. 


2 . 6 iLffect of Domain Structure 

In deriving (2.5) it has been assumed that the static 

magnetic field is large enough to break down the domain structure 

and the specimen can be regarded as a single domain magnetized 

to saturation so that the anisotropy effects might be neglected 

In tne case of Weaker magnetizing fields there will be a 

multidomain structure in the specimen as a consequence of which 

additional peaks will appear m the absorption spectrum. 

6l 

Nagamiya who calculated the effect of domain structure on 

FhR m single crystals with tetragonal symmetry has explained 

6 2 

the additional peaks observed experimentally by Bickford . 


2 . 7 Phenomenological Approach t.o Damping 

The divergence of dynamic susceptibility /f x to ® 
at resonance (60 = (x} Q ) means only that attenuation (damping) 
processes have not been taken into consideration. It has been 
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pointed out in the preceding chapter that a part of the width 
of the ferromagnetic resonance line arises out of the inter- 
action of the uniform mode with other modes such as phonons, 
magnons and conduction electrons of the system. The effect 
of these interactions can be taken care of formally m the 
phenomenological theory by adding to the right hand side of the 
eqn (2.1 ) the relaxation term wmch describes damping 

The first phenomenological equation was suggested by 
Landau -Lifshitz^ 7 m 1935 and has the following form 

= - /(lax's) - M x {7* x 1?) (2.13) 

dt M * 

where ?) is the relaxation (damping) parameter 

A slightly different way of writing this equation was 
48 

suggested by Gilbert whose equation reads 

-# = - Y a <i?.b) ♦ $ («*#) 


The basis of writing t is form for the relaxation term 
m ( 2 . 13 ) and (2.14) is that it represents the torque which 

— 5 

characterizes the damping of the magnetization vector M to 
B along a spiral path. 

Another phenomenological equation was introduced by 
Blooh to describe nuclear magnetic resonance and was employed 
by Bloambergen to explain ferromagnetic relaxation. The 
Bio ch-Blo ember gen^ phenomenological equations read 


dM 


dt 


x.y _ 


dM z 


- /(«? * 3) X ,y 

- /(l? xt) z - 


M 


z 


M 


X 



t 


1 


(2.15) 


dt 
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where C ^ and are relaxation parameters r 1S the 

saturation magnetization This form of the relaxation term 

characterizes the exponential decay of the transverse components 

of magnetization to B with a characteristic time constant 

For weak radiation fields when we put M = M , the 

z s 

Landau-Lif schitz-Gilbert and Bloch-Bloembergen equations are 

equivalent The shape of the absorption line m all these 

cases is jl, orentzian 

64 

Yager et al studied the ferromagnetic resonance m 
single crystals (spherical samples) of nickel ferrite (NiO FegO^) 
at room temperature at 24 KffiC/sec The value of the g-factor 
obtained by them was 2.1 9 and the first order magnetic anisotropy 

1 1 O 

constant was found to be -6*27 10 erg/cnr as compared to 

4 

-6.2*10 + 10$ observed from hysteresis loop measurements by 

65 

Galt et al . They employed Landau-Lif shitz damping and were 

able to fit the resonance curve quite satisfactorily. 

A similar FMR study on small spherical samples of nickel 

66 

ferrite was made by Healy over a wide temperature range 
between -195°C to 588°C. The g-factor and first and second 
order anisotropy constants were obtained from the study of 
the variation of the lmewidth with temperature. The g-value 
was found to be essentially constant over the entire temperature 
range. The anisotropy constants and showed different 

sign behaviour over the temperature range, Healy was, however, 
able to fit his resonance curves by either type of damping - 
Landau-Lif shitz or Bloch type. His results showed that the 
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lmewidth decreases with, the .inc rease m temperature which 
suggests that the spin-lattice and dipolar interactions do 
not contribute materially to the lmewidth. The measurements 
on spin-lattice relaxation times made by Damon^ have indicated 
that the spin-lattice interaction contributes little to the 
lmewidth. The conclusion drawn by Healy that the spin-lattice 
relaxation does not contribute to the line broadening is not 
totally correct as the relaxation time occurs m the phenomeno- 
logical equation as a gross effect and cannot be separated in a 
simple manner. 

68 

Young et al who investigated FMR in nickel and 
supermalloy m two different conditions of magnetization, namely, 
( 1 ) the applied magnetic field is parallel to the surfaoe of the 
sample and (2) it is perpendicular to the surface, found that the 
shape of the curve was correctly predicted only m the case of 
parallel magnetization and m the case of perpendicular magneti- 
zation there were marked discrepancies (in the low field region) 
such as the appearance of a second minimum m the frequency 
shift curve. This they tried to explain on the basis that in 
the porpendicular case they were in the region of unsaturated 
magnetization and the theory is consequently unable to explain 
the phenomenon. They found that though both the Bloch and 
Landau-Lif shitz damping terms lead to Lorentz shape of the 
curve, the former was more satisfactory than the latter. 

According to the usual theory of FMR, the ratio of 
absorptions in the two cases of parallel and perpendicular 
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magnetizations is independent of the numerical value of the 

2 7 

damping constant hut its dependence on the parameter 

68 U) ° 

is different m the two cases Their results differed from 
the predicted ratio m the most favourable case (parallel 
magnetization case) by a factor of about 1.5> an< i they were led 
to interpret that neither type of damping term was appropriate. 
Ament and Rado^ gave an electromagnetic theory of 
"spin wave resonance" experiment and satisfactorily explained 
the exchange shift as well as the shape of the absorption line. 

They used Landau-Lif shitz type of damping m their calculation. 

70 

Rado and rfeertman' , who studied spin wave resonance m 
single and polycrystals of nickel iron alloy having extremely 
low magnetic anisotropy at 3 and 4 KkC/sec, showed that the 
resonance curves could be interpreted on the basis of the 
theory of Ament and Rado ^ by using an effective exchange 
stiffness constant A of (3*3 + 0.5) x 10 erg/cm and a 
spectroscopic g factor of 2.06 + 0.01 They found, however, 
that the agreement (qualitative) could be destroyed if a 
Landau-Lif shitz or Bloch damping term were included. 
Experimental data taken at liquid nitrogen temperature showed 
an increase m the resonance lmewidth and a decrease in the 
resonance field A part of the discrepancy was attributed to 
the anomalous effects caused by the electronic mean free path. 
Kore recently Hirst and I range 71 have shown that the discre- 
pancies encountered by Rado and tfeertman in trying to fit 
their nitrogen temperature data cannot be attributed to the 
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non-local conductivity (anomalous skin depth effect). 

The destroying of the agreement by using either type 
of damping ( Landau -Lif shit z or Bloch) is not very surprising, 
Aado and vfeertman^ studied the effect of exchange on the 
lmewidth, and their experimental conditions were such that 
the exchange term was predominant and the relaxation term 
was made negligibly small by suitably choosing a specimen of 
zero anisotropy, of large conductivity and by working at low 
field Therefore, addition of any damping term will spoil 
the agreement with the experiment, as the relaxation processes 
were not operating. 

. .31 

The ferromagnetic resonance lmewidth m a metal is 
believed to be due to ( 1 ) magnetic relaxation damping whose 
microscopic origin is not exactly known. it is described by 
including Landau -Lif shit z or Bloch-Bloembergen type of damping 
terms in the equation of motion for magnetization and (2) 
exchange conductivity broadening which arises due to the fact 
that the exchange stiffness opposes the tendency of the dynamic 
magnetization to vary m the skin depth It depends upon the 
exchange stiffness constant, conductivity and magnetic aniso- 
tropy that might be present It has already been pointed out 
xn the preceding chapter that the lmewidth due to exchange 
conductivity broadening varies as (where Co* is the 

frequency) and that due to magnetic relaxation damping as Lx). 
The measurements on FMR linewidths in nickel and iron have 
been made by many workers 33 ' 35 . These authors have used 
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Landau-Lif shitz type of damping Their results are not all 
m agreement with each other. Salamon 3 ^ found that the value 
of the g-f actor m the case of nickel suddenly increases as 
he goes from below T c to above T Q (T being the Curie tempera- 
ture) , while Bhagat et al 33 observed that the increase in the 
g-f actor was rather gradual. The reasons for this discrepancy 
are not known It is possible, however, that the spin-orbit 
coupling effects might be very important m the case of nickel 
The band structure effects of the d-electrons and the influehce 

of the conduction electrons on them should also be taken into 

72 

consideration. Very recently Kambersky has studied FMR m 
iron whiskers and has reported lmewidths at the highest 
frequency ( 73 GHZ). The lmewidths found by Bhagat et al 52 ’ 53 
at the corresponding temperatures and frequencies (upto 35 GHZ) 
are larger than those obtained by Kambersky. 

2 . 8 Discussion 

On the basis of the interpretation of the results of 
measurements on FMR using phenomenological theory, one oan say 
that the th§ory explains the shape of the resonance line and 
also the anisotropy satisfactorily. The lmewidth and its 
temperature dependence is not, however, explained by the theory 
adequately. The reason is probably that as absorption is a 
gross effect and the relaxation termsused also represent an 
overall effect due to all the damping mechanisms present in 
the system, en d so the relaxation parameter can be chosen to 
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fit the resonance curve nicely, but the details of the line 
cannot be reproduced by such phenomenological parameter whose 
microscopic origin is not exactly understood 

We, therefore, conclude by saying that the phenomenological 
theory is not adequate for a quantitative description of FMR 
lmewidths . 
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CHAPTER - III 

GREEN FUNCTION THEORY OF FMR AND MAGNQN -PHONON COUFLING 

3* 1 Green Function Formalism 

Solxd state theory problems are inherently many-body 
problems The usual quantum mechanical perturbation theory 
methods are inappropriate to treat these many-body problems. 
However, in recent years, there has been a considerable interest 
in trying to treat these problems by the methods of quantum 
field theory. One of these metnods is the method of double 
time temperature dependent Green functions, which has been 
applied very extensively and successfully to a wide range of 
problems of solid state physics such as ferromagnetism, super- 
conductivity and lattice dynamics, etc. The method is non- 
perturbative. The poles of the Fourier transform of the Green 
function yield information regarding the elementary excitations 
of the solid and their decays. Fost of the quantities of 
immediate physical interest can be obtained in terms of these 
Green functions. 

,tfe can define different kinds of Green functions the 

double time causal Greenfunction defined in terms of the average 

value of the T-product of operators, or the retarded and 

advanced Green functions. Throughout the present work, we 

shall use retarded Green functions There exist several 

reviews^ and texts'^ on Green functions. W« would closely 

73 

follow an excellent review by Zubarev . The retarded Green 
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function is defined as follows 


W tft * = B(t' )3)= <^A(t), B(t')J> (3.1) 


where a and 3 are Heisenberg operators 


A, B = e' 


A, B 


-lHt 


H is the system Hamiltonian 


0(t) = 1 


for t > 0 


= 0 otherwise 

(a system of units m which ii = 1 is used) 

The bracket \ > ~] denotes the commutator and X - N denotes 
the grand ensemble average, namely, for any operator 0^ 

°p = T r ( eX H O p )/T r e-# H 


Here T r denotes the trace of the quantity that follows 
= (jI^T) , i ^ being the Boltzman constant and T is the 
temperature of the system The operators A(t) and B(t^) 


satisfy the equation of motion 


C a, hJ 


(3-2) 


Differentiating the Green function (3*1) with respect 


to time we obtain 


id<j( t-1 
dt 


= «((')]■) * ^(t),H3 , B(t ' p. 

(3*3) 


The Green function <^^(^)» ^ J > xn ec l n * (3»3) 
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general involves Green functions of higher order than the 
original Green function B(i/ )^> , except, of course, 

for the cases of non-interacting systems where exact solution 
can be obtained. One is, therefore, forced to linearize the 
eqn (3*3/ by suitable decoupling approximations 

The time correlation functions which are of physical 
interest are defined as 

F 3A (t,t ^ = B(t/)>(3.4) 

In the statistical equilibrium, these functions F^ A and 
depend, as do the Green functions, only on the difference 
t-t^. Equal time correlation functions are called auto- 
correlation functions and they give the average values of 
the products of operators One can express the correlation 
functions in terms of the spectral density function J ( (*> ) as 


\ 

-O 

1 

-p 

= au>j(co) e - 1 ^(t-t / ') 

(3.5a) 

-CB 



7°° J ,, \ ft 60 -i A)(t-t ; 

= J du)J(w) e' e 

(3.5t>) 


-CD 


Since G( t ■ 

-t ) is the difference of F AB and Fg A , it 



follows that 


G(oo + iL) - g(CJ -i£) = -i(e^ -1) J(&>) (3.6) 

E. — ^ 0. 

Here G{6J ) is the Fourier transform of G(t-t ) defined 
as follows 

a(a ,) ' r *') oft-/) d(t- t') (3.7) 

A __ 

— 00 
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The expression (3.6) relating the Green-function and. 
tiie spectral density function is important Cnee the retarded 
Green function is known, the spectral density function can he 
obtained and then one can get various correlation functions at 
all times and at all temperatures 


3 . 2 Linear Response Theory 

Kubo^ has given an elegant formulation cf the problem 
of the response of a quantum mechanical system to an external 
perturbation. The theory assumes that the response of the 
system to the external influence is linear. The higher or the 
nonlinear terms are neglected m the above analysis, but they 
are important m the theory of phenomena which arise due to 
nonlinearity. 

Let us consider a system with time -independent 
Hamiltonian H q . A perturbation H^ nt which depends on time 
explicitly has been switched on at t = - r c, when the system 
was in thermal equilibrium. Now the system is described by the 

Hamiltonian 


H 


H 


H 


int 


( 3 . 8 ) 


The ensemble average of any operator A is given by 

<A > = »[/(*)*] (3 ' 9 

nUera f(t) is the density matrix representing the ensemble 
and satisfies the following equation 
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= H H ' ^<*>3 - [>„ +H* nt , f( t)] (3.10) 

Let J (t) = J q + ^^(t) such that 

^/( t ) = 0 for t = - cd and /(-oo) = f = — — 

7 '° Tr e-/ 811 

(3.11) 

Restricting to the linear term m the external perturbation 
and making use of the initial conditions (3*11 )i the change m 
the density matrix can be obtained as 

^f(t) = -x /* dt / (3-12) 


Now the average value of the operator A is given by 


A = ,rr L {A - 1 A J 

.OO J _ y 

= A o + 1 J < A (t), « int > dt (3.13) 

-OO ' 

t ' / / 

Considering = F(v ) 0 , where F(t ) is the force 

coupled to the operator part 0 of the dynamical system. Then 

j > *h v - ±H ,y 

H t = F(t ) e 0 e ° = F(t ) 0(t ) (3-14) 


Then the effect <&A = A - A q will be given by 

<AA = J G(t-t ) F(t ) dt (3.15) 

Thus the retarded Green function is essentially a response 


function. 
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3 3 Dynamic Mag netic Susceptibility 

In studying ferromagnetic resonance one has to find the 
changes m the average magnetization components 
(i = x, y) caused toy the radiation field B x (t), as they yield 
us dynamical magnetic susceptibility 


^ m (t) = X B(t) (3.16) 

whion determines the transfer of energy from the radio frequency 
field ^hG spin system Remembering that 


Am x * ^ S |(‘) (3.17) 

--4 

nere f refers to the lattice site. One can see that the 
problem reduces to evaluating the average changes m the spin 
operators ^S^(t), which is independent of f' because of 
translational invariance 

The Hamiltonian of a ferromagnetic system m the 
prosence of a radiation field can toe written m the form 


H 


H 


o 


+ H 


t 

int 


(3.18) 


Here H ^ describes the interaction of the spin system and 

the radiation field, and H includes the interaction with 

0 

external static magnetic field and all other interactions 
such as exchange interactions and dipolar interactions, etc 
Assuming that the constant magnetic field is directed 
along the z-axis and the oscillating field lies m the x-y 
plane and is uniform throughout the sample, we can write for H 
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S^t) s?(t), (, 

J > f 


x, y) 


(3-19) 


It is obvious that in the case of a plane polarized, radio 
frequency field and m the linear approximation, we can write 
for ^ the following expression 


X 


= -2 


.2 y r 

'f 


b ^ 4 , s r 4 


% 


( 3 . 20 ) 


l h 

ij / 

Thus to evaluate bhe dynamical magnetic susceptibility we must 
find the Fourier transform of the Green function . 


3.4 Hamiltonian of the System An Applicati on of Green 
Function Method 

Consider a ferromagnet m a constant magnetic field B 
in z -direction. The Hamiltonian 56 can be written m the form 


H 


where H z 


H 


ex 


and 


H + H + H, 


z 


ex 




(1 = x,y,z) 


( 3 . 21 ) 

( 3 . 21 a) 

(3.21b) 

(3.21c) 


? — * . 

Hero all the symools have their usual meaning J(f 1 -f 2 ) is 

. th 

the exchange integral which is positive, Dy denotes 
component of the magnetisation and ^ is the demagnetisation 

factor . 


The three terms in the Hamiltonian (3.21) respectively 
denote Zeeman energy, Heisenberg exchange energy and classical 
demagnetizing energy. The demagnetizing energy (3.21c) arises 
due to the finiteness of the sample dimensions and, therefore, 



42 


takes into account the demagnetization effects of the surface 
of the finite sample. 

The demagnetizing energy (3.21c) can be written m the 
quantum-mechanical form using the relation 


M. 


■ /“b i- r s f- 


(3.22) 


It will be convenient to change over from the spin 
operators to Pauli operators 

S~. = , S = i(c£* _ Cj>), Sj> = 1 -2ng*, ng> = Cp C^(3 23) 

where C^' and Cp are annihilation and creation operators 
for spin deviation at the f th lattice sitre. These operators 
obey following commutation relations 

CpCf = o , Cp Cg' - = o , f £ f' 


(3.24) 


C - ctt> + C-t C — 1 ss 1 , C& C = Cp Cp = o 

f f iff* ii XI 

Using (3,22) and (3.23) the Hamiltonian ( 3 . 21) can be written 
in the following form 

H = -| N J(o) - B + £ 2 B + 2J(o) j % n f -> 


-2 


z op c f . - 2 z, 

"f^ft, 12 *1 *2 ?lf 2 




) n f° 

X 1 £ 


+ * N » a /4 - N, z: , "r + “f, 

,2 

, 4-2 /V „ N ) 7 . r Cjrt 0 ~> + ci* c5r> *1 

+ 2N V ? L°f 1 f 2 f ! f 2 J 


N_ 


n. 


2 

Ab 

■fa m i — 

N 


< N x * V 4? 0f 

12 


2 


(3.25) 
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*-i -J 

where J(k) = 23, J(£) e lk ‘ , J(o) = o 

f 

xs the Fourier transform of J(f) 

From equations (3*20) and (3 23) we obtain 

2 XX - c f> + 4.°r °f*> 

+ cgog, o f> ] 


(3.26) 


To determine /j , we should evaluate the Green 
/l/ XX. 

functions appearing in (3*26). Using eqn (3.1), we define 


the Green function G^.-* through 


G sf" = ^ C g > C f^ 


(3.27) 


gf 

and write down its equation of motion 
dG-s„> 


i ' - s S# S(t-*' ) * 1 2 ^ B - 2 /4 N * + “ (o) J Ggf " 

-2 - 4 s < N * + V J G p f,+ 2 ^l? J(8 ' p,) 


- <*, ♦ N y) ] «»,<*. 4? -»£(«(**> 

*-4 °S- (3 ' 28) 


where *-= 1 - 2 4 ng> « tbe ” latlVe -8 netiaat1011 

Thie equation show. that higher order Green function, hays 

a i 'TVf-a •wi’i'fc© down thcs© h.ig'iior 
made their appearance as usual. T 

of the original ones , we 


order Green functions m terms 



kh 


c f> 

“ 11 G pf 

°F>> 

= ^^pf’ 


use the following decoupling procedure^ 

< n r c r °J*> - 

c ? ^ ^ ‘^ n g'> °f* ^ 

» tos V = <<#. c r > 

With the above decoupling approximations (3*29) the 
equation of motion of the Green function G ■ ^ reduces to 

if - S(t-t ) ♦ <-&■ ("x-’V) L.V 


(3 29) 


dG 




x y‘ 


* B + 2 J(o) -2/4 N z fjv 

< N x + N y)] G pf (3-30) 


-Z ^ r 2J (s -p) r - IT + "y 1 J u pf 

Tp 1 


2 

ft-B 

N 


In the same approximation the equation of motion for 

F „ can be written in the form 

g f 

i Zlf' _ (n -n ) Z fly - f 2 A 3 B + *~\ 2J(o) 

1 dt n v x y ; 4=£* pf L / B u 

2 

-2/4 N z J] ^ + L 2J(? - S) - (WJ 1 *' 

(3*31) 

Introducing the following 


and 


(3*32) 


fa") ^ 


w(t-t' ) 


d ^ 
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Using eqns (3-32), ( 3 . 30) and ( 3 . 31) we obtain the 
following set of equations for the Fourier transforms of G 
and F functions 


< ^ " L %) “ S k . j0 <r ^4 (N x - N y ) F£(Aj) = 


cT 
27 T 


wher 


£•£,0^/4 ( N x ~ N y) G k(^) + ( ^ + = o 

o S jg -- 2/Ug (B -tr^U B N z ) + <T(2J(o) - 2J (lc) ' 

+ ^/4 < N x + »y) ^,o 


(3-33) 


(3.34) 


We shall solve eqn (3«33) separately for the case 
lc ^ o and the case k = o. 

Case 1 k ^ o 

If lc £ 0 we get from eqn ( 3 • 33 ) 


( - 6*) Gg(CO) = i.e., G^(U)) = 27r (£j-f R ) ’ F ^ = ° 

and ( 3 • 35 ) 


Sj^ * 2 


/fe (B -/yU B N z ) + 2*"{j(o) - J(k)j 


Here is the energy of an elementary excitation 

m a ferromagnetics "spin wave energy" Because of the factor 

df , the elementary excitation energy depends on the 

temperature . 

f/fe can use ‘^ ) to evaluate (T* • We know that 

JC 


n s= 


<V> - <°r c r> 


Now using the spectral density relation 
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n 


±y r 

k -ext 


oo 


e 1 _1 


%;( (-o+ i£) _ Gg(to-x£jjdtu (3. 


and oqns (3.35) and (306) we obtain 

1 r 7 cr" 


n s 




N - e 


-1 


From this we got 


Ph 


» 5 2j Ooth ( C_I£ ) 

13 N -jr’ v 2 ' 


(307) 


Tho quantity £.^>, as pointed out earlier, determines 
tho ©J ©montary oxoitation spectrum m the system and cannot 
be negative (if we assume that can take on negative 

values, wo are lod to a contradiction* n < o). The solution 
(305) is thus possible only if ° for a11 As 

tho second term in £.£ is non-negative, we can get the 
solution considered here only in sufficiently high fields 


B - <r 


/ U B N z ^ 0 


(308) 


which moans that the external static magnetic field must 
"overcomo" the demagnetizing influence of the boundaries, 
i.e,, the ferromagnet should behave as a single domain. 
Case 2 T? a o 

For He » o, we get from oqn (303) 

cST~ 

(W~ e o ) G 0 (CO) (N x - N y ) F o (o>) = jyr 


(309) 


( N x “ Ny) * 0 (U) + £ o> 


= o 


36) 
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c 

0 

2 

* 

Xb b 

+ <<74 («x 

+ N - 2N ) 
y z ' 

(3.40) 

From 

the 

above 

equations , 

we get 


0 o (w) 


a" 

+ c 

0 




27 r 

x ,2 2 

^ “^0 



'o< w > 




- N ) 

y ; 

(3 41) 


'2 rr 

2 2 




i 


where ^ 0 =/A^ = 2^ Qb h tfyig(N x -N y ) j B +dyU B (N y -N B )?J 

This is called the resonant frequency for the uniform precession 
mode and was first obtained by Tyablikov 56 . 

Using the relation 


< A - B > u =-< A+ - B+ ^ 


(3-42) 


and eqns ( 3 . 26 ) and (3,4l ) we can write down the following 
expressions for the real and imaginary components of the 
dynamical magnetic susceptibility ^ 

1 


/XX 


Xo ( A a/ A i F 


1 - 


03 * 

~~~2 

U)o 


and 


where 


% 


n 

XX 


la - 


«7 4 


(3.43) 


(3.44) 




In writing eqns (3*43) and (3*44) we have used the 
following identity 


It 


x + iC 


» p ~ + 17 r£ (x) 

A 
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Here P denotes the principal value 

/' 

The imaginary part ~Y of the dynamic susceptibility 

/wcx 

character lies the absorption of energy by the spin system from 
the radiation field In the approximation considered here, 
this absorption occurs at a frequency of the r.f field equal 
to W Q - 

3 

Kittel was the first to obtain the value of the 
resonance frequency using the equation of motion method He 
obtained the following expression 

h, - *m, [i B + <v n z > m H 3 * <v n z > m }]* {%h5) 

The resonant frequency expression ( 3 . 4l ) obtained by 
Green function method is a generalization of Kittol's formula 
(3.45) and is valid at all temperatures. 

From eqn (3.41 ) it is clear that for an unbounded 
medium, i.e., for a ferromagnet of infinite dimensions 
reduces to 


60 o = Z f*B B SinCS = ° 

The resonance frequency, therefore, evidently depends 

on the shape of the specimen It has also been established 

, .. 60,62 

experimentally 

Discussion 

We notice that m the present approximation the 
IfwWvgnetxo resonance line turns out to be delta-function 
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shaped This is connected with the fact that in the present 
approximation the Green functions have poles on the real axis 
In the higher approximation (spin wave - spin wave processes) 
the Green functions would have a cut along the real axis and, 
therefore, one would get damping It should be pointed out 
that the coupling of magnons with other excitations of the 
system such as phonons and conduction electrons would lead to 
lines of finite width 

In the next section we shall consider nn .non-phonon 
coupling. The coupling of the magnons vith conduction electrons 
will be considered m Chapter IV 

We have considered here only the demagnetizing effects 
of the surface of the finite sample. The volume demagnetizing 
fields have not been included. Actually the dipolar interaction 
can give rise to volume demagnetization which produces 
elliptical precession instead of circular precession. The 
ellipticifcy is maximum for magnons, i.e., when the angle 

between the spin wave vector and the applied static field is 
90° Haas 76 used the dipolar Hamiltonian and calculated the 
renormalized spin wave energies The results obtained by 
him for uniform precession mode using Callen decoupling 
procedure are found to be identical with those of the spin 
wave theory and extend them through the entire temperature 

range . 

6 Magnon -Phonon Coupling 

It has already been pointed out in Chapter I that m 
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ferromagnetic insulators the dominant relaxation mechanisms 
are due to magnon-magnon and magnon-phonon interactions. The 
magnon-magnon interactions arise out of the exchange and dipolar 
interactions. Various workers have investigated these processes 
extensively in garnets. The magnon-phonon processes originate 

1 4 

from the magnetic anisotropy, magnetostriction and dipolar 
interactions 1 ^. The observed magnetostriction can be accounted 
for by a free energy density of the form 


f = V* 2 e yy * f 9 «) + B 2(^“xy + / 3 /V + /V e zx ) 

(3.46) 


where B 1 and B_ are constants to be determined from the 
experiment o( . and / denote the direction cosines 

of the magnetization and e's denote the elastic strain 
components The strain components can be obtained by 
differentiating the displacement with respect to appropriate 
spatial coordinates. The displacement u ± (of the atom or 
ion at 1 th lattice site) can be expressed m terms of phonon 
creation and annihilation operators bjr, and b-.as follows 


V — — » 

ziy^r - iK.r, 

- _ (3.47) 

where v = (^)* £ j S. [4^ + b -l > + J 

Here V is the volume and ^ is the density of the material, 

£ denotes an orthogonal set of polarisation vectors, C. 

' ' n „ o a (c = C is the longitudinal 

fdehote the phonon velocities (G^ °l 

A 

| r _c is the transverse velocity). 

JreJLocity and C 2 


(3.47) 



Using (3*46) and (3*47) the following expression for 
the Hamiltonian in the second quantized form can be obtained 
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H = Fk k L E(K) ^ a H + E * (K) ak ' ^ 1 (3 -“ 8) 

k=k' +K 

where i5(K) = constant (K)^, describes the coupling of magnons 

with phonons, a* and a^ denote respectively the creation 

and annihilation operators for magnons of wave vector k The 

terms in (3.48) represent two magnon - one phonon processes. 

The first term says that a magnon of wave vector k and a 

phonon of wave vector K annihilate to create a magnon of 

wave vector k. The second term denotes the inverse process, 

i.e , a magnon of wave vector k is annihilated to create 

anothor magnon and a phonon The energy and momenta are 

conserved m these processes There may be other processes 

such as one magnon - one phonon which have been ignored since 

78 7 ® 

they contribute insignificantly' in YIG. Kasuya et al 
have shown that one magnon - two phonon relaxation processes 
are also unimportant in YIG. The magnon phonon relaxation 
processes are not very important in the case of metals 


3.7 Hamiltonian for the Coupled Ma g non-Phonon System 

Considering a system of magnons interacting with phonons 
the Hamiltonian for such a system can be written m the form 
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where H, 


m 


H 


and 


H. 


ni-p 


-2T ^k^ a k ®k + 

^ C k ^ b K b K + 

27 L E(K) ^ ^ + E ‘< K > \ 7 

- > «*£ t Ja. ' 


<^ - p y 

< 

k 

k=k T +K 


\ = 


N, 


= ^ b : 


'K 


O 

Here 65^ = gyktg B + Dk is the spin wave energy and = k, 

is the phonon energy D is the spin wave coupling constant 

The various terms H , H and H _ m the Hamiltonian 

m p m-p 

(3*49) denote respectively magnon, phonon said magnon-phonon 
Hamiltonians. a+, and b+, bg have already been defined, 

n^ and Ng represent the occupation number operators for 
magnons and phonons respectively E(K) which has already been 
defined describes the strength of the coupling of these two 
excitations . 


3.8 Magnon and Phonon Green Functions 

We define single magnon and phonon Green functions as 


follows 


<Cei- ^ q (*)- 

<CV 4) - < (0) > 


(3.50) 


tTa jng eqn. (3.3). we write down the equations of motion for 

*i 

Fourier transforms 
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F ; +K , K>q ( .) +¥ ^ (K) Wi ,m 

“ d (aj - e ^ )G - (W) w„(^) 

vhere «e have introduced the following mired Oreen functions 


q+K,K,q 


q-K,K,q 


(60) 


(CO) 


<Vk*k’ <> 

CO 

< a q + K b K> a q > 


(3.52) 


T k + q,k,q(^ ) ‘ <C a k+q & k ’ b q 

h) 

We construct the equations of motion for these mixed 
Green functions and get 

(U - *W + 8 K ) F q + K,K,q - 2 '(‘K a q+K+ k *£. »; ^ 

+ Z Ok) ^ q+K _ k k* b k , a^>(3-53a) 

- E *(K)r <^< +K a k a q+K , 

^"H-K * V F q-K,K, q = S E ( k ) ^ a q . K+ k \ b k’ & q ^ 

+E ( K )J^J<^_ a k a k+K a q~K* a q^> ( 3>53b ) 

'CJ 

o^L ^ E ‘ (k] «V«\ V 

,(^-6^. +q +6/ k ) T k +q}k>q (^) = "^ E ( k 1 )<^k+q a k-k 1 b k , b+j^> 

1 'OJ 

+ 27 E( k ik^a k+ b+ , 

1 l 3* 53c) 

+ ^ i E " (k 1 ) ^< a k + q. ki ^ b ; 

*k + q b q 


&V eV 
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These equations (3*53) do not contain an inhomogeneous 
part as K^o , i.e., E(k)=o for K=o. 

We notice that eqns. (3*53) contain still higher order 
Green functions for which we should, m turn, construct 
equations of motion and continue the process further We do 
not wish to attempt this here and decouple the higher order 
Green functions as follows 


< ^ a q+K+k 

b + i 
K 




G m 

q+K+k,q 

(") = o 


<C a q + K-k 1 

b + 1 
K 

b k' 


\ ? Kk G “q <■ 


(3- 54a) 

<C a k+K a k 

a q+K’ 


“q+K ^qk 

G qq 



<C a q-K b K 

b k 

, a + 

q 

** 

( 1 +S K ) ^Kk G qq ( ^ 1 



^ a k a k+K 

a 

q 

-K’ 


a q-K^k,q 

„ g“ (W) 
-K qq 


(3.54b) 

<^ a q-K-k 

b K 

b k’ 

a + \s - 

q C 

< b K V G 

m tu)) - 

q-K-k, q 

o 


^^q a k 


k i 

’ b *> 
*> 

“ <$k+q a k- 

*,) «% ’ 



! ^. a k+q+k 1 

n + 
°-k 

b k 

k i 


” sfk+q+k 1 


b q» ■ ° 

(3.54c) 

^> a k a k+q 

" k 1 

b k 

k i 

• b «> 

" “k ^qk 1 

! 

G P tu}) 

qq 




and a k+q b^ , b+^, ^ “k +q ^qk 1 G qq ( - U)) 

.There n^. = <a* and Ng. = b K ^> 

In the above equations the correlations between two 
Ration and two annihilation operators have been neglected 
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With these approximations (3*54), eqns (3»5l) and 

(3.53) Yield 


G m (JO) = 
qq v 1 


1 


zTr\co- -<W)] 


and G*L(A/) = 

qq 271 ^ - 


L - M J 

q q 


("0 


(3.55a) 


(3 55b) 


where 


M“(AJ) =Z j®(K)j 


N IC ~ n q +K 


1+N K + n q-K 


U) ~ <^q+K + & K ^ ~ ^q-K - k 


- «"> ■ M'5 


K 


(3.56a) 


(3.56b) 


The poles of Green functions G qq W and G qq twjrespectlV ® ly 


are given by 


and 


uj- cj q -*£(«) 
W - E q - Mq(<^ ) 


Using the identity 
1 


(3.57a) 

(3.57b) 


C — * o 


x + 10 


P - 

X + 


ifl< 


f(x) 


We write down for the real and imaginary parts of 
u/) and M p (<^) respectively 

Qi — 


£(«) - He m“(W) = pZ|e(K)| [J-u**x+“k * "-"^q-K- G K 


UN K * °q-K 


(3.58a) 
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+ (,+S K + Vk^W-V^ 1 ] (3 58b) 

H q ( "> = Re M%(fr) = P|E(q)| 2 ^_ T (3.59a) 

k lc+q k 

and 

/q (W) = Ira ■?("> ■ 1 £ M 2 z ( fi k -a k «)5(^-6/ Jt4 ,*W k )(J.39b) 


where i indicates that the principal value of the corresponding 
integral should be taken 

The functions OJ) and y^(£v/) play the role of 

m „ 

damping for magnon and phonon excitations. P (Oj) and I l y (CJ ) 

Q. Q 

denote the corresponding energy shirts for them 

Using eqns. (3.6), (3*55), (3«58) and (3.5 9) we can 


obtain for magnon and phonon distribution functions 


m R^ 

1 r - 1 ) 

~irr“ -Jr i 2 . 


and N 


r [v-V'rfi*')] 

1 j m /v 1 ^ ^ dw 

7T - \u. E ^(-yjj 2 +[j%UJ)] 


( 3 . 60 a) 


(3 60 b) 


If ), y^(4J) and and R^( 6i ) are 

Slowly varying functions of &J and if and 

;Ve smaller than r”‘(c/) and R q ( U ) respectively then 
eqns (3«6o) assume the following form 



(3.6la) 
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and 


N 


cJ 

f* 8 , 

' q - 1 


V w 

where and £ axe determined by 

q h 


and 


v 

to q 

rJ 

e 


/V 


£ - M^(£ ) 

q q v q' 


( 3 . 61 b) 


- - <(^ q ) - 


= o 


(3 62 a) 
( 3 . 62 b) 


We notice that n and N are Bo se-Emstem 

q q 

distribution functions for magnons and phonons which are "smeared 
out" because of magnon-phonon interaction The widths of the 
"smearing out" regions are of the order of the magnitude of 

rV 


the 


damping ^q ( ^ q ) and Vq(^q) respectively. 

The elementary magnon and phonon excitations are given 

79,80 




ft/* /y/ 

P, 


by 

tJ 


K 


l z r 

\" n q+K 

1+S K + “q-K 1 


(Mf 4> q+x + c K ) 



and 


ft/ 

£ = L 

q q 


♦ !E(q)|PZ' 


k (£ q" 0 k+q 


(3.63a) 


(3.63b) 


We observe that the magnon and phonon energies are 
renormalized as a result of the magnon-phonon interaction 
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CHAPTER - IV 

MAG-NETI C REL AXATION IN METALS "s-d EXCHANGE 11 INTERACTION 

k . 1 Spin Waves m Metals 

In the Heisenberg model of a ferromagnet besides the 

s P iri aligned state constituting the ground state, there exist 

spin wave states close to the ground state For metals where 

the electrons are assumed to be itinerant, it was not clear 

for quite sometime whether spin waves could exist in them 

The existence of spin waves m both ferromagnetic metals and 

ft 1 

insulators has been established experimentally . Herring 
82 waves 

and Kittel who proposed that spin/did exist m metals showed 

that the magnetization of a metal responded to a magnetic field 

with a sinusoidal spatial dependence in such a way as to 

indicate the possibility of spin wave like excitations 

In an itinerant electron model, the electrons are 

divided between 'up' and ’down* bands. These ’up' and 'down* 

spin bands are displaced relative to each other due to an 

exchange interaction which acts just like an internal 
8 3 

magnetic field . Consequently, the wave vector for the 
electron in the highest filled state will be different for 
•up' and 'down' spins. An electron in the 'down' band can be 
excited to ’up* band without any change in the wave vector 
x>f the electron, provided it is given sufficient energy, 
yhile an electron from ’up’ band can only make a transition 
jbo ’down’ if its wave vector changes by at least -^k^, the 
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horizontal spacing between the bands 

These transitions are described as single particle 
excitations as we reverse the spin of a particular electron. 

Using all the electrons, a collective excitation can be obtained. 
This is done for the 'down' to 'up' transition m which, as 
stated earlier, the wave vector change may tend to zero and 
it is this excitation which is assumed to be a spin wave m the 
itinerant electron model. The wave vector of this spin wave 
will be equal to q, the change in wave vector of the electrons 
as they go from one band to another 

We discuss spin waves m an electron gas which is 
assumed to have a net polarization. 

8 4 

The Hamiltonian of the system can be written in the form 

H -S t k * Z . V(ql (£ +q> ^ C*._ w .C kV C k ,,(4.1) 

where C + . and C denote creation and annihilation operators 
k o k « 

for Bloch waves of wave vector k and spin ^(o.) ^- s ^h.® 

Fourier transform of the Coulomb potential 
We define a quantity S k (q) 


S k< q > = C k+q T °kl 


( 4 . 2 ) 


Which describes an excitation m which an electron with down 
spin in the state k is removed into a state k + q with spin 
up. The equation of motion of this quantity can be written as 


dS"(q) 

dt 


js k (q), H*] 


(*• 3 ) 



The commutator of 


60 


S k( q ) with the Hamiltonian (4.1) can be 

evaluated to yield 

dS'(q) 


«r* 

Ck + k 2 ), C ^<f ~ ^k+q+kg r- C k 1 -kg £ C k 1 1 ( 4 *4) 

¥e apply random phase approximation (rta) tc the right 
hand side of eqn, (4.4) and retain only such terms m the four 
operator products as contain pairs of operators C + , C and 

Id P K ff 

replace them by n^^, where is the mean value of the 

number operator ^ in ground state. The equation 

of motion (4.4) now reduces to 
dS“(q) -7 

1 dt - (C k- £ k tq ) + £ ' , ( k 1 )[/ k+ki+q ,-* k+ki Js;(q) 

+ (^k* - *k + <H V < k 1> s ; + k( q > < 4 ' 4 > 

Introducing the renormalized energies ( J j 1 through 

id P 


£ k- e k+q ) 


=k 


(<i) + TL v(k )[ 


k 1 1C 2 


C + C + 

k 1+ k gj/ k+q^ 


= & k ” ^ V ( k l) Vkd 

K 1 1 

^kt = £ k " 23 v ( k -|) 


(4-5) 


and assuming that the ground state is ferromagnetic so that 
the numbers of the up and down electrons are unequal and so 
are the energies 6^^ and ^kj,’ we kaVe 

( i ft + ^k+qt' W k i) S k( q ^ + ^“k+q'fr ~”k-P 23 v ( k l) S k+k/ q ^ = 

( 4 . 6 ) 


o 



6 1 


If the states of the 


system are assumed to be of the form 


l l 'b~ A k s k<i) i°> 

wo-th energy CO, then the HP A equation for the coefficients 
Aj^ can be written as follows 



It xs obvious that such states, if they exist, will 
resemble spin waves since they reverse one unit of spin 


As mentioned earlier, V(k 1 ) is the Fourier transform 
of the Coulomb interaction. However, if correlation effects 
are taken into account the screened Coulomb interaction will 
be more appropriate. Replacing its Fourier transform by a 
constant V (V "J o) we obtain 


1 + v z ( *W - l 

k k*k+q^ “ ~ ^ 

The relations (4.5) now reduce to 

" £ k - N * V 

and <O k t * £ k “ N t V 


(4.8) 


(^* 9 ) 


where N^, and denote total number of electrons with 

down and up spin respectively. We assume that )> Ity and 
consequently <✓ LO^. 

We notice that the seoond terms in the right hand side 
the eqns . (4,5) and (4.9) denote the exchange energy of 
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the electron It appears that the Coulomb interaction is 
responsible for displacing the spin bands for up and down 
electrons. However, some authors are of the opinion that 
this displacement of bands is caused by the exchange inter- 
action 

For vanishingly small value of q the eqn (4 8) 
takes the following form 

1 + vr 2^ ] = o (4 10) 

^k^F (N-V " N^)V - 6J 

^k+q T /C p 

Using the eqns. (4.9) and (4.10) we obtain 


1 


V 


\ + 1r< y v 

s k 7 vv 


v ( N t . up - ! k+q ♦ E k + W 


(4.11) 


or 

1 5 = 


^ ' Nj ' ) E k^ E F + V + 


23 


E k > W 


V(N t -N^)* 


< > 


We now have 




7— 1 7% V k) °k + iU-^k) 2 ^k * •] < 4 - 12 ) 

(N^ -Nj,) U 

If £ = £ the first term of (4.12) vanishes and 


wo get 


i <( q. V k ) 2 ^ k > 


(4.13) 


frhere the average 4^- >is taken over the region between the 
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two Fermi surfaces 

This establishes that there are low lying excitations 
of the spin wave type . Thus we have the same spectrum as 
spin waves for the localized spin model. It should be pointed 
out that m this calculation V(q) has been replaced by a 
constant V greater than zero, which means that V(r) m 
real space is a delta function This is not unreasonable m 
that the interactions are strongly screened by conduction 
electrons. tfhile this particular type of sinqplif ication is 
not essential to bring out the general nature of the excitations, 
it do os show tnat some type of localization is brought into 
the model 

Rajagopal et al*' 0 have studied single particle and 
colleot d vo excitations m an itinerant electron model More 
recently, Sokoloff 85 has also investigated polar spin waves 
in ferromagnetic metals. 


4.2 "a-d Exchange 11 Interaction 

The spin wave system can also relax through the 
interaction with the conduction electrons This interaction 
can be brought about by the "s-d exchange” interaction. This 
is regarded as the most predominant interaction m ferromagnetic 


transition metals. A number of workers 41 ' 44 have studied 
this interaction in detail In the present work, this 
interaction has been studied to explain relaxation in these 
metals. It ha. already been pointed out in Chapter I that 
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the exchange interaction between 3 d-electrons is mediated, 
by conduction s electrons Actually, owing to tne "s-d 
exchange" interaction between conduction electrons and the 

j. 

localized 3 d spin at x lattice site, the conrucfcion electron 

spin system is polarized and consequently the spin polarization 

*kliX 

of conduction electrons appears around the 1 site Usually 

th 

most of the induced spin polarization is localized at a 
site and contributes to self energy. A part of the induced 
spin polarization may extend to other lattice sites 3 and 
induce an effective exchange field on the Jd electron spin 
localized at j th site through "s-d exchange" interaction 

4.3 Hamiltonian 

The Hamiltonian for the 3d-localized ferromagnetic 
electrons exchange coupled with conduction s electrons is 
of the iorm 


« - H z + H C + H ex + H s-d 

(4,14) 

where H z =- & n 

(4 . l4a) 

H *2 b k C k <r C kr5' 
c kd" 

(4. 14b) 

H - - X 

ox n >m 

(4.l4c) 

and H. d - - S Vd 

5 -u mn 

(4. l4d) 


Tha various symbols have their usual meaning. The 

Side of eqn. (4.14) denote 


three terms on the right hand 
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respectively Zeeman, conduction electron, exchange and 
»s-d exchange" energies. 3h terms of the creation and annihilation 
operators c i Z( f anc * c k/( tf* = + -j^the "s-d exchange" Hamiltonian 
(4 . l4d) can be -written m the form 


H 


s-d 


1 

N 

nkk 


) e 


i(k'-k).R n , 


ic. c, S + c, c . f S_ 
Ik- k + n k+ k - n 




C, r -C. C. 

k - k+ k 




l4d ) 


where S— = + i S^, J .(kk ) is the s-d exchange 

ii xi •*« ns "~ct 

integral given by 


J s-d^ kk/ ) = Ijy lj r ) 0*( T ') M r ) d3p d3r ' 

whoro $(r) is the atomic function for d-electron and 
is the Bloch electron wave function. 

86 

Introducing Hostem-Primakoff transformation and 

using spin wave approximation, i.e , S n = b^ , s n = b n 

where b + and b create and destroy respectively spin 
n n 

deviation at the lattice site n, the Hamiltonian (4.14) 
can be expressed as follows 


H 


\ ■ ( f>* l? q 1 

£(k -k+q) [o; + V _ b q * <£'_ c k+ bqp S 5 

S(k -k+q'-q) o fc -_ -4 + V.] b q V (4 ' 15) 


where b + and 


q 

if or spin-waves 


b are creation and annihilation operators 

q 

of wave vector q, and 
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8 J 8-d (q ‘ 1 } + 2 ^ B B ' 



2 

q_ 

2m 




The 3rd and 4th terms represent the coupling of magnons 
Vvith conduction electrons. The 3rd term denotes the 
scattering of an s electron from one Bloch state to the other 
including a spin flip, and simultaneous creation and destruction 
o X a spin wave The 4th term represents s electron scattering 

0 

without spin flip and spin wave scattering from state q to q 


4 . 4 Magnon Green Functions 

We define a single magnon Green function 


G (t) B 

/in ' ' 


qq 


<**(*>. *;,(•)> * -•*<*> <[y %1> (4 - 16) 

Hero, as has been pointed earlier, ^ denotes averaging 

over a grand canonical ensemble at a given temperature, 

©(t) = 1 If t>o, and zero otherwise Here all the 
operators are in Heisenberg representation. 

The time derivative of ^qq^ ( * ) 13 given 

dG qq. 1 (t) n c , , (!S)K J ,(k,k-Q.) X 

i __! . b(t) a qq , + ^ 8 - a 

< c iU- °k t - ’0 + ^ ^ Js - d(kkl } * 

< ( c £_ °k 1 _ "°k+ °k, b ktq-K, ’ ^ 4 ’ 7 ^ 

Here we notice appearance of higher order Green 

„ 0 +Vlo following decoupling approximations 

functions. We use the following 


to decouple them 



67 


C°k- 

C k 1 - b k+q-k 1 ’ 

V 

”k- ^ldc 1 

G , 





(4 18 ) 

C lc 1+ b k+q-k 1 ’ b q 1 ^> 

N# 

”k+ ^kk 1 


where 

5 k- = <T C k- C k~ > ’ 


- <<♦ 

c k+^> 


Fourier transformation of ( 4 . 17 ) with equation ( 4 . 18 ) 

yields 


[w -V 5 Sj B . a (^.k) (\. 


2 rr 


-(f ) i ^ J s-d(' [ ' k ^) F 1 c. q ,k, qi ( w ) 


where we have substituted 



(t) 


_ 1 s . t , V l Wt 

N j dco G (<a>) e 

qq., 


(4.20) 


“ d F k-q,k,q 1 (W) = 'Cfk-q- V’ b q, ^ 

To obtain G n _ (w) from oqn ( 4 . 19 ), we should write 

H.H -j 

down tho equation of motion for ■ F i s >q,k>q,j ^ ^ ^ 

derivative is given by 


*Vq,k 

1 dt 



fS + - L~ ) F (t) “( iP - x 

v & k k~q' k-q.q.q^ N 

(Z d s-d( k - k - p) C°k-q- V*- V b q,> 

L -o 


- s d s .d( F ' k - q) 


p+ c k+ 


1 





s-d 



o 


+ Z J 8 .d( klI, K 0 k-q- V 

PP-J 


b p + q-k’ b q^3 

. C k+ b p J b p o +p-k+q’ b q 

b P 1 b t+p 1 -p ,b q 1 ^>j 

(4.21) 
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The so equations contain severs! higher order nixed Green 
functions which are decoupled in the following way 

C k+ b p+q-k’ ~ S pk G qq ^ ft) 

^°P C - °k + b J- b p 0+ p-keq-^ K. S P k- q ,k,q W ^ 

<°k-q- °k-p- b p ' ^ )) « Vq- §qp G qq^> 

“ d <°k-q- V b p, b k +Pl -p' b q,) i S Pi r k.,,k,q 1 vt'*y i *-2=) 

•“* J, 

n p = 4 . b ^ > is the average number 


where n - 

P 

of spin waves , 

Taking the Fourier transform of eqh, (4.21) and using 
eqns (4.18) and (4,22) we get 


F. 


/2S,4 J s-d^ k,k ' q ^(”k-a- " "kJ 

k-q,k,q^ ) = ^ N ) 7 77 71 2^* c- . 


Here 


s ~d 


U>- 87 + 6" + | j , « qq 1 

k k-q N s-d p 

P 

Vd( k " q * k ~ q ) + J s -d^ k,k ^ 


From eqns. (4.19) and (4,23) we obtain 

Sqq 1 / 2 7T 


G nn ( W ) = 


qq 


^ 0 - c q - A- f^) 


(4.24) 


where we have put 


S ZJ J s.d( k ' k ) K- - 

, 2S V J a-d< k,k-q ) ("k-q.- ~ 

/ 4 (Cd) . T ig 

' W - + E k-q + N J e-d Z -V 


(4,25a) 

(4.25b) 
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we get 


Separating the real and lmaglnary parts of » , (J) 
get J H 


k U-S, + +8“ j X f 
k k-q N s-d p *p 


(4.26a) 


)^(Ch) = Im /q(W ) = -^ J Ld(k . k . q) (5 k . q _-H k+ ) a 

5(^-6; IS' + | ? s _ d X n ) (4.26b) 

P 

Here n p denotes the superposition of elementary Bose 

excitations with resonant intensity and damping \/(U)) t 

' Qi 


and is given by 


i r 

n = - J 

p /T J 

1 -oo 


[U-8 q - ^^(60)j 2 + ^( W )-j 


(4.27) 


We know that the response of a ferromagnet in a 
static magnetic field and subjected to a radiation field 
will be given by the magnetic susceptibility tensor 
defined as follows 


ty M (t )> 


(4.28) 


where M(t) is the magnetization at the time t. For the 


present case of a ferromagnetic metal it is obvious that 


M * ®d j^B ^d + S s (4.2 

where S, and S are the total spin operators for d and 
d s 

s-electrons respectively. Using (4- ..28) and (4.29) we have 


(^.29) 
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^ V M y^ - A‘b [*4 < s d' ^ + S s e d ^ sj, s* ^ 

+ * 0 ® a « S I> s^ +g J <3 y, s ^| (4-30) 

The magnetization of s-electrono ■, a Dra( ,n 

electrons 13 small compared with the 

ma gne tization of d—electronsi mi, _ * , 

trons The latter is always appreciably 

larger and specially at small values of the ratio J g A 

being the Fermi energy), the magnetization of the s- F 

eloc trons can, therefore, be neglected. For the present case, 

the resonant properties of the system will be described by the 

Green function /£ S y . S Y 

d 1 d * 

In terms of the spin wave creation and annihilation 
operators Y becomes 

/lyy 

Xyy = g2 ¥ [<tv b o + <C?o’ b o>/ ( 4 ' 31 ) 

From the above expression we notice that at resonance the 
spin waves of k = o should be excited in the framework 
of the spin wave theory. It is called uniform precession 
mode ferromagnetic resonance, ^ is obtained using 
eqn. (4.24) for q = o as follows 


V .2 2 SN | 60 ~ ^o 60 + C o *) 

A" Y b ^ L + ^ f£ 0) + (w+ e o ) 2 + /^)i 

. /£ SN ( 

+ 1 *{(u>-l 0 )\f 0 (G» - (U.t 0 )\yl^) 


r*f 

Ia) + c 


g 2 SN 


(^. 32 ) 


e 0 + A + Re f 0 («) 


jEJere 
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and 


x< 


4) 


) = 


2S 

N 


2 J *-d( k .k) ("k--”**) S(C0-f£+Cfc+|j s _ a 2jn ) 


The use of <§.A + , B + ^ = - <^A,B ^ * y has been made an 

obtaining eqn (4.32). From this equation, it is clear that 
the line is of the Lorentzian tyre However, it is known 
experimentally that at resonance the spin waves excited possess 
^ ^ 5 being the skin depth) of the order of reciprocal 

skin depth. The width, therefore, should be given by 


A<A = 


2 K^) 

h 


(4.33) 


where ) ls given by eqn (4.26b) and has the 

significance of the lmewidth. It is obvious from the expression 
for ^(W) that for finite non-vani shing lmewidth the 
following energy relation should be satisfied 


co - si 


V 


'k-q 


Z7 S 


- 
n 


-d ^ 


- o 


(4.34) 
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The last term m eqn (4.34) is an additional term and is 

43 

absent from the corresponding equations of Vonsovskn. et al 

The above energy conservation relation suggests that 

in the relaxation process a bunch of magnons are participating 

43 

VonsovskiL et al J have drawn attention to the fact 
that the relaxation process which results in spin flip of a 
conduction electron by emission or absorption of spin waves 
becomes ineffective for low energy spin waves due to the large 
exchange gap between bands of s electrons with different 
brientations and energy conservation law. However, 
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Heinrich et al have recently suggested that non-stationary 
states of s spins make possible, m principle, a non-resonant 
absorption at a very low spin wave energy It is quite possible 
that the thermal magnons denoted by eqn (4 2J) might be 
associated with such non-stationary states of s-electrons 

It should be recalled here that the ferromagnetic 
3d-electrons have been assumed to be localized and a Heisenberg 
type of exchange has been used to describe them The 3d- 
electrons in ferromagnetic metals, specially Ni, are not 
localized but are known to be itinerant Therefore, an 
itinerant electron model should be a realistic approach to 
the situation existing m these metals 

It will be interesting to mention here that recently 
88 

Hubbard and Jam have reported a calculation of the 
generalized spin susceptibility in a correlated narrow 
energy band model which is based on the narrow energy band 
theory given by Hubbard®^ According to this theory, a 
partly filled narrow energy band (s-band) is considered and 
it is assumed that the electron electron Coulomb repilsion 
acts only between electrons of opposite spin on the same 
atom and has the magnitude I. They have discussed two 
special cases of their susceptibility function, namely, 

(l) weak interaction (band) limit, i.e., I where & 

is the band width and (2) strong interaction (atomic) limit, 

i.e*, I 

In the wealc interaction limit their result reduces to 
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yo 

random phase result obtained by Izuyama et al and in the 
opposite limit a Curie law is obtained. Their theory is also 
not likely to explain lmewidths m ferromagnetic metals. 

It has already been pointed out that the spin wave 
energy is negligibly small compared to the energy required to 
flip the spin of the conduction electron. If the spm wave 

3 4 

damping is to take place, as many as 10 or 10 spin waves 
would be required to participate m the relaxation process. 
This seems to be an impossibility However, such processes, 
if they operate m these metals, should be reflected m spin 
wave resonance experiments. The superconducting ferromagnets 
may, perhaps, be the best materials to look for such processes 
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CHAPTER - V 

CONCLUSION AND DISCUSSION 

In the present work the interaction of the ferromagnetic 
3d electrons m a ferromagnetic transition metal (Ni, Fe) 
interacting with the conduction s electrons by means of "s-d 
exchange" interaction has been considered to explain ferro- 
magnetic relaxation. Using Green function formalism the 
expressions for the dynamic magnetic susceptibility and line 
width due to "s-d exchange" interaction have been obtained. It 
has been suggested that a bunch of magnons participate m the 
relaxation process resulting m the conduction electron spin 
flip. It should be noticed here that an exchange interaction 
of the Heisenberg type has been assumed between the magnetic 
and conduction electrons This is not justified as the 
Heisenberg model assumes that the electrons are localized on 
the atoms at the lattice sites The electrons m transition 
metals are not localized but are known to be itinerant (band) 
electrons. This is because of the fact that the 3d shell in 
these metals is on the outer side of the ion and, therefore, 
the electrons can move from one atom to the other forming a 
narrow band of states. A band model for magnetic electrons 
would, therefore, be a more realistic approach to the situation 
existing x» ferromagnetic metals and the effect of conduction 
, electrons should be properly taken into account. 
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The spin wave excitations in. an itinerant electron 
model have been discussed. It is found that the spin wave 
energies in the long wavelength limit are proportional to 
k (k being the wave vector of the spin wave) and thus have 
the same spectrum as the spin waves on the localized electron 
model. We thus have a probable explanation for the use of the 
localized electron model to a ferromagnetic transition metal 

The renormalized spin wave and phonon energies resulting 
from the coupling of the spin waves with phonon* have been 
obtained. Xt is found that the Bose -Sms tern distributions 
for spin wave and phonon occupation numbers are "smeared 
out" because of roagnon-phonon coupling 

It should be mentioned here that the "s-d exchange" 
interaction can be used to investigate many other interesting 
properties such as conduction electron spin polarization, 
anomalous electrical resistance of transition metals and a 
very interesting phenomenon, i.e., the coexistence of 
ferromagnetism and superconductivity in some superconducting 
ferromagnets , Xt should be noticed that "s~d" or "s-f" 
(interaction between the conduction electrons and the 4f 
unfilled shell electrons) will act unfavourably to the 
simultaneous realization of the ferromagnetism and super- 
conductivity. If the process of the spin wave damping 
results in spin flip of the conduction electrons then the 
i»s-d exchange" will oppose the tendency of simultaneous 
‘PfQ.Gurnehc© of ferromagnetism and superconductivity , as the 
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absorption or emission of spin waves by, say, one of the 
electrons of the Cooper pair resulting m its spin flip 
will break the Cooper pairs and, therefore, will destroy 
superconductivity. 

The theory discussed here is unable to explain lmewidth 
and its temperature dependence m ferromagnetic metals The 
reason for this is that it has not been possible to take the 
exact microscopic origin of the relaxation mechanism into 
consideration. In order to find the microscopic nature of the 
relaxation mechanism, more refined experiments on FMR 111 
ferromagnetic transition metals must be performed, which can 
give some clue to the mechanisms m these metals. 
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